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Comment on "Shadowability of Statistical Aver- 
ages in Chaotic Systems" 

Lai et al. 1] investigate systems with a stable 
periodic orbit and a coexisting chaotic saddle. They 
claim that, by adding white Gaussian noise (GN), 
averages change with an algebraic scaling law above a 
certain noise threshold and argue that this leads to a 
breakdown of shadowability of averages. Here, we show 
that (i) shadowability is not well defined and even if it 
were, no breakdown occurs. We clarify misconceptions 
on (ii) thresholds for GN. We further point out (iii) 
misconceptions on the effect of noise on averages. 
Finally, we show that (iv) the lack of a proper threshold 
can lead to any (meaningless) scaling exponent. 

(i) Shadowing deals with macroscopic error propaga- 
tion of noise bounded by a small value (e.g. computer 
accuracy of 10 -16 ), comparing a noisy and a 'true' 
trajectory 0. Since the authors use GN, one cannot 
properly speak of shadowing. However, even for bounded 
noise, their claims are still doubtful. In a recent work, 
the error of an average was found to be amplified by a 
factor of up to 10 12 , although the trajectory never exits 
the attractor This renders a reliable computation 
truly unfeasible. In 1], though, the effect on the averages 
is only of the same order as the variation of the noise 
level D. Furthermore, the averages in Figs, la and 2a 
in do change even less below D = 1CP 2 ' 5 . Therefore, 
shadowability is not compromised at all. 

(ii) The authors of |l( mention a threshold for GN 
above which the periodic orbit and the chaotic saddle 
would become connected. Yet, such a threshold does 
not exist, since the mean first exit time (t) is given 
by Kramers' law (r) ~ exp(A^-), where D is the noise 
variance and AU is either the potential [J] or, for 
nonequilibrium and chaotic systems, the quasipotential 
difference 0. Thus, (r) varies with D, yielding a 
different average for every noise level. 

(iii) Generally, averages depend on the noise for all 
noise levels, implying that no threshold can exist, even 
with only one metastable state. For the linear map 
x n+ \ = ax n + b + with a fixed point = and 


white GN one gets (x 2 ) 


D. This is pictured 


in Fig. la, fitting perfectly the data. Although the 
average appears to be constant for low noise, it depends, 
in fact, on the noise for all D. The same applies also to 
nonlinear systems (cf. fit in Fig. lb) and bounded noise. 

(iv) Because the value of the threshold is arbitrary, any 
scaling can be achieved-just by tuning D c , as fittings of 
the form of Eq. (1) of [l| are very sensitive to the value 
of D c . To verify this, we show in Fig. lb the logistic map 
x n+1 = ax n (l - x n ) + D £ n with (£ n ,fm) = <W as in 0. 


The putative threshold of D c = 10 5 is marked by an 
arrow. The graph is evidently not constant there. With 
D c w 7 • 1CT 6 (first arrow), a scaling AG ~ (D - D c ) a 
results, with AG of 1]. This yields a w 2.1 (Fig. lc), 
in clear contrast to a rj 1 reported in |l|. Therefore, no 
reliable (i. e. any) exponent can be obtained, since D c 
is not well defined from the outset. 
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FIG. 1: (a) (x 2 ) for the linear map vs. D with a = 0.7 
and b = 0.75 (circles) and analytical fit (x 2 ) — x 2 + 1 _ 1 aii D 
(full line), (b) (sin x) for the logistic map with a — 3.8008 
and a polynomial fit for D < 1.4 • 10 _J (full line) (c) Scaling 
AG ~ (D-D c ) a with D c « 7-10" 6 (circles) and least squares 
fit (full line). All averages are for 10 8 iterations. 


As to the reply, the authors claim to have "a theoreti- 
cal justification for the existence of a threshold" through 
D c = y/ A$/ lnx -1 , with A<I> the quasipotential (see 
(ii)) and \ the probability resolution. Contrary to Ref. 
[5] of the reply, where \ drops out since only propor- 
tionalities under parameter variation are considered, D c 
in depends for fixed parameters on \ & n d hence a is 
not uniquely defined. Applied to (iv) with D c = 1CP 5 of 
and A$ w 7.8 ■ lCT 10 (not shown) yields X = 0.0004. 
But a finer resolution x = 10~ 7 gives D c = 7 ■ 10~ 6 , 
with a « 2.1 (cf. (iv)), whereas \ = 10 -2 results in 
D c = 1.3 • 10~ 5 and a « 0.21 (not shown). Again, this 
invalidates any proper scaling. 
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